3D-printed ultra-small Brownian viscometers

Measuring viscosity in volumes smaller than a microliter is a challenging endeavor. A new type of microscopic viscometers is presented to assess the viscosity of Newtonian liquids. Micron-sized flexible polymer cantilevers are created by two-photon polymerization direct laser writing. Because of the low stiffness and high elasticity of the polymer material the microcantilevers exhibit pronounced Brownian motion when submerged in a liquid medium. By imaging the cantilever’s spherically shaped end, these fluctuations can be tracked with high accuracy. The hydrodynamic resistance of the microviscometer is determined by fitting the power spectral density of the measured fluctuations with a theoretical frequency dependence. Validation measurements in water-glycerol mixtures with known viscosities reveal excellent linearity of the hydrodynamic resistance to viscosity, allowing for a simple linear calibration. The stand-alone viscometer structures have a characteristic size of a few tens of microns and only require a very basic external instrumentation in the form of microscopic imaging at moderate framerates (~ 100 fps). Thus, our results point to a practical and simple to use ultra-low volume viscometer that can be integrated into lab-on-a-chip devices.

The second large group of micro-rheological techniques employs vibrating cantilevers immersed in liquid.AFM (atomic force microscopy) cantilevers of ca.100-400 µm length are commonly used.The viscosity and the density of the medium are determined by evaluating the resonant frequency (typically in the range of 10 3 -10 5 Hz) and the quality factor of the oscillators 25 .Cantilever vibrations are stimulated either actively, via external actuation [26][27][28][29][30] , or passively, via thermal fluctuations 25,31,32 .Different MEMS-based cantilever-type viscosity sensors were also reported 33 .
In this work, we introduce a new type of microscopic viscometer whose working principle is based on the passive Brownian motion of a low-stiffness flexible microcantilever.The microcantilever is created by two-photon polymerization microfabrication, with a support microstructure anchored to a microscope coverglass.A micronsized spherical bead is formed on the cantilever's free end, allowing for highly accurate imaging-based tracking of its thermal fluctuations in a liquid medium (Fig. 1).The viscosity of the surrounding fluid is then obtained using the power spectral density of the measured cantilever fluctuations.The micro-viscometers were tested and validated in water-glycerol mixtures of known viscosity.

Viscometer fabrication
OrmoComp photoresist and OrmoDev developer were purchased from Micro Resist Technology.SylgardTM 184 from Dow Corning GmbH, Germany was used to prepare PDMS layers.High-purity distilled water was used to prepare water-glycerol mixtures of known viscosity.The viscosity of the mixtures was calculated using an empirical formula 34 .
The flexible structures were fabricated by TPP-DLW (two-photon polymerization direct laser writing) of OrmoComp using a 785 nm, 100 MHz femtosecond laser (C-Fiber 780, Menlo Systems, Germany).A 40 × oilimmersion objective (NA 1.2, C PlanApo, Zeiss, Germany) was applied to focus the laser beam into a photoresist droplet, which was placed on top of a 0.17 mm thick glass coverslip.A closed channel with two inlets made of 4 mm thick PDMS block was formed over the droplet.An acousto-optic modulator (APE 130,101, APE GmbH, Germany) adjusted the laser intensity.The laser focus was scanned by moving the sample with a piezo stage (P-563.3CLPIMars, Physik Instrumente, Germany) in three directions.The microstructure design is shown in Fig. 1.Four identical cantilevers were placed on a single support.The nanowire beams were drawn vertically as a single line, and their stiffness was adjusted by setting the scanning speed and the laser power.Optimal conditions, ensuring both flexibility and self-support, were obtained at 10 μm/s scan speed and 5 mW laser power.After the polymerization process, the uncured photoresist was washed away by filling the channel few times with OrmoDev developer and ethanol.Just before developer-ethanol exchange, the coverslip was irradiated with the mercury lamp of a fluorescence microscope (HBO 50) to finish the polymerization and stabilize the microstructures 35 .Finally, the microfluidic channel was filled with distilled water.
Microstructures with designed bead radius, R d , of 2 µm, 4 µm, and 6 µm were tested.The length of the cantilever was 20 µm in all cases.The sphere centre was positioned at 15 µm above the coverslip surface.The typical diameter of the used cantilevers was 0.35 µm.

Material properties of photopolymer microstructures
The elastic properties of TPP-DLW-fabricated microstructures have already been investigated.Displacement recovery measurements of stretched coil-shaped micro-springs made from SCR500 and PMMA photoresins were reported [35][36][37] .Most importantly, the bulk and microscopic shear moduli of the photopolymer material were found to differ significantly, with microscopic values being orders of magnitude lower in some cases.The Young's modulus of photo-polymers in micron-sized objects were measured by deflection tests using nano-indentation equipment 38 .The elastic moduli mentioned thus far only describe the (near-static) slow motion of photopolymer microstructures.When operating the microstructures in a dynamic regime involving Brownian fluctuations, the viscoelasticity of the photopolymer must be considered.The micro-viscometers investigated in this study are made of OrmoComp, a bio-compatible photoresist material with exceptional flexibility 39 .Our previous work investigated the viscoelastic properties of OrmoComp nanowire cantilevers by measuring the recovery motion of structures deflected by an optical tweezer setup 40,41 .It was shown that the cantilever mechanics can be described by a three-parameter viscoelastic model.Furthermore, an analytical solution for the power spectral density of freely fluctuating viscoelastic cantilever/microsphere structures immersed in Newtonian liquids was obtained 42 .After being refined, this result serves as a theoretical basis for the presented viscometers.

Fluctuation measurements
The viscometers submerged to different water/glycerol mixtures were viewed in brightfield under a microscope equipped with a camera operated at 663 fps.Brownian fluctuations were captured by taking 30-s videos.The viscometer fluctuations were assessed by image tracking the bead position.We used a custom MATLAB algorithm that employed a direct search scheme to determine the spatial shift of the bead on each recorded image with respect to the first recorded image.Using normalized cross correlation (normxcorr2 MATLAB function), the algorithm first found the spatial shift with 1 pixel resolution.In the next step the algorithm probed half pixel shifts in each direction added over the already known 1 pixel resolution shifts.For this, we spatially shifted the image with fast Fourier transform, using Fourier Shift Theorem, and calculated the mean absolute pixelwise difference between the shifted images and the reference image.The algorithm then selected the smallest mean difference value among the nine half pixel shift combinations (for directions: ± X and ± Y, and 0 shift) and repeated the search probing quarter pixel shifts.This step was repeated eight times to achieve a resolution of 2 −8 ≈ 0.004 pixels, corresponding to 0.34 nm.Only the one-dimensional data along the Y axis (see Fig. 1) were used for further analysis.

Data analysis
We begin to derive the theoretical background of the current viscometer operation by defining the fundamental postulates that mark the model framework.The first assumption concerns the viscometer inertia.The effective Young's modulus of polymer nanowires is four orders of magnitude lower than that of silicon, which is commonly used in AFM and/or MEMS cantilever applications.The stiffness of nanowires employed in our viscometers is extremely low, on the order of 0.001 pN/nm 41 .At the same time, because of significant dissipation caused by the surrounding liquid and the cantilever material's internal viscoelasticity, low frequency motion dominates the system dynamics.The cantilever and the attached bead motion is strongly overdamped with no resonances observed in the studied low-frequency range.Due to these facts, all inertial forces and the related effects are neglected in our study.Owing to the cylindrical cantilever shape, the two possible oscillation directions X and Y are assumed to be independent and identical.The second assumption concerns the viscometer geometry.The external viscous damping force of the surrounding liquid is assumed to act primarily on the attached bead, representing the thin nanowire limit.This latter assumption will be refined later by using an empirical formula to correct the model for the finite cantilever thickness.
Figure 1d depicts the mechanical scheme of the micro-bead attached to the viscoelastic cantilever; only a single oscillation direction is considered.The nanowire viscoelasticity is characterized by the upper (red) arm, which evolves from the Kelvin form of a three-parameter standard linear solid 40,41 with two spring constants (stiffnesses) k 1 , k 2 , and a single internal damping coefficient δ.The lower (blue) arm containing the hydrodynamic resistance γ represents the Newtonian liquid that damps the micro-bead fluctuations.Based on this model, and taking the previous assumptions into account, the power spectrum of the observed Brownian fluctuations obeys a double Lorentzian distribution 42 : The two corner frequencies f c1 , f c2 , as well as the auxiliary K A parameter are non-trivial functions of the three viscoelastic parameters (k 1 , k 2 , δ) and the hydrodynamic resistance γ (see the Supplementary Materials; for more details we refer to 42 ).In the frequency region above the two corner frequencies, the power spectrum simplifies significantly: The criterion to reach pure 1/f 2 dependence is to stay at least one order of magnitude above the two corner frequencies.This frequency region is identical to the high-frequency power spectrum limit of a bead trapped in a harmonic potential 15 and reflects free diffusive Brownian motion 43 .We emphasize that the expression of Eq. 2 is independent of the cantilever viscoelastic properties; the hydrodynamic resistance γ is the only structure-related parameter that affects the fluctuations.As γ is proportional to the solution viscosity η, the amplitude of the power spectral density in the 1/f 2 regime can be used for viscosity measurements.A single calibration parameter, ξ cal , is introduced to relate the solution viscosity to the measured hydrodynamic resistance: In practice, a calibration measurement in a liquid with known viscosity (η ref ) can determine ξ cal for a specific viscometer structure: (1) www.nature.com/scientificreports/Two additional effects are considered when fitting the experimental power spectral density data in the theoretical 1/f 2 region.First, due to finite sampling frequency the experimental spectra are aliased near the Nyquist frequency 15 .Second, video capture introduces a blur effect that needs to be corrected 44,45 .Fitting the experimental power spectral density data by the corrected theoretical function is required for further data analysis.To reduce experimental noise in the power spectral density, the Welch method with n = 32 non-overlapping time windows is used.The experimental curves are analyzed using least-square fitting in the 20-90 Hz spectral range, with theoretical values used as weights (for details we refer to Ref. 46).Finally, the obtained power spectra are corrected for the bias introduced by least-square fitting when combined with spectral averaging 46 .

Results
Figure 2a depicts typical fluctuations of microcantilever structures immersed in water.The observed pattern has a minor asymmetry, as shown in Figure S1 of the Supplementary Materials, which depicts the position probability distributions for the X and Y directions.The difference in distribution widths (around 13%) is probably due to the inherent asymmetry of the polymerization laser beam, which is then transferred to the cantilever.The time course of the bead fluctuations along the Y axis (plotted in Fig. 2b) is used for data analysis.The probability density of the bead center position is fitted by a Gaussian distribution (Fig. 2c) as expected for Boltzmann statistics.The deviation from the Gaussian distribution at the tails (see also Figure S1) is most likely caused by a minor drift in the camera image.However, a low frequency drift has no effect on the fluctuation power spectral densities in the frequency range used for data analysis (see below).Using the equipartition theorem, the distribution variance is used to calculate an effective trapping potential stiffness 44 .The obtained value, k eff = 0.8 10 −6 Nm −1 is of the same order as the equilibrium stiffness reported previously for similar viscoelastic cantilevers 41 .The Supplementary Materials include a video of a fluctuating 4 µm radius bead.
Figure 3a depicts the measured power spectral densities for three different bead radii (R d = 2 µm, 4 µm, and 6 µm).The power spectra are relatively noisy in the low-frequency region, making the data fitting with double-Lorentzian functions over the entire frequency range difficult (see Eq. 1).However, determining the parameters  www.nature.com/scientificreports/ of the double-Lorentzian fits is not the goal of this work.Rather, we focus on the frequency region (above 20 Hz) where theoretical power spectra decay with the second power of the frequency (see Eq. 2 and the dashed line in Fig. 3a).The higher the bead radius, the lower the power spectral density in this region.Due to aliasing, the measured power spectra level off near the Nyquist frequency (332 Hz). Figure 3b depicts power spectra measured with a R d = 4 µm bead in various water-glycerol mixtures.A higher viscosity (higher glycerol content) decreases the power spectral density in the 1/f 2 domain.This effect can be used to measure viscosity, as shown below.
The power spectral data were fitted in the frequency range 20-90 Hz, as indicated by the grey regions in Fig. 3.The hydrodynamic resistance γ obtained from the fit for different solution viscosities and different bead radii is plotted in Fig. 4. As expected, the measured γ values exhibit a linear dependence on the solution viscosity and can therefore be accurately fitted by linear curves that pass through the origin.The experimental points are well aligned along the fitted curves, which paves the way for viscosity measurement applications and allows for a simple calibration procedure.The slope of the fitted curves equals the calibration constant reciprocal, 1/ξ cal (see Eq. 3).In principle, a single measurement made in a solution of known viscosity is sufficient to obtain the value of ξ cal for a specific microstructure type.Larger bead radii yield smaller calibration constants.In the Supplementary Materials we provide a more detailed analysis of the relationship between the calibration constant and the bead radius.

Discussion
Flexible stand-alone microstructures composed of a nanowire cantilever with a spherically shaped free end were introduced to measure viscosity of sub-microliter volume samples.The power spectral density of the microbead Brownian fluctuations decays with the second power of the frequency in a region spanning more than one order of magnitude, which is consistent with the proposed theoretical model.As expected, good linearity was found between the hydrodynamic resistance obtained by fitting the power spectra and the solution viscosity.The microstructures were calibrated using the known viscosity of water-glycerol mixtures.
The primary benefit of evaluating power spectral densities in the 1/f 2 domain is that the data in this region are unaffected by microstructure mechanical properties (see Eq. 2).This is especially significant when we consider that the mechanical properties of photopolymer microstructures can change in different environments due to the penetration of solute molecules into the polymer pores, as previously demonstrated for aqueous glucose solutions 41 .Fortunately, the 1/f 2 region can be easily reached using the low-stiffness flexible cantilevers presented in this work.Surface adsorption of analytes should not interfere with measurements if the structure dimensions (bead diameter) do not change dramatically.It was found sufficient to evaluate data in the 20-90 Hz range, where the limiting power spectral density defined by Eq. 2 is well above the noise level even at the highest viscosity studied.The 30-s videos recorded during our measurements represent a good balance between the required data storage capacity and measurement time and the good quality fits to the power spectral density data.
The Brownian micro-viscometers presented in our work have a different operation regime than other viscometer types that rely on thermal fluctuations.In comparison to microbeads trapped in optical tweezers, a simpler optical setup consisting of a microscope equipped with a fast camera is all that is required to measure viscosity.Furthermore, the current microstructures are prepared at well-defined positions on the substrate, e.g., in microfluidic channels, so there is no need to disperse them into small-volume samples; the structures are easy to localize under the microscope.In comparison to AFM and MEMS cantilevers, we highlight the small dimensions of the polymer structures.The stand-alone viscometers tested here have a characteristic size of ~ 30 µm, which can be further reduced.There are no additional (bulky) supporting parts required to operate the viscometers.In principle, sample volumes well below one microliter are possible and sufficient for viscosity measurements.Some technical issues must be resolved to turn the current structures into a commercial product.One possible development path points to disposable chips with pre-fabricated microcantilever viscometer structures.

Figure 1 .
Figure 1.The viscometer structure.(a) Brightfield image of microstructures with beads of 2 µm, 4 µm, and 6 µm radius.(b) A screenshot from the video (see the Supplementary Materials) shows the 4 µm radius bead.(c) Schematic view of the solid support with four vertical cantilevers.A spherical bead is attached to each nanowire beam.(d) The mechanical model of the viscometer with the viscoelastic beam and the microbead immersed in a Newtonian liquid.

Figure 2 .
Figure 2. Brownian fluctuations of the microbead attached to polymer nanowire.(a) Typical record of the R d = 4 µm bead center position in water as tracked in the X-Y plane.(b) The time course of the bead center fluctuations along the Y coordinate during the 30 s data acquisition period.(c) The probability density of the bead center position fitted by a Gaussian distribution (σ is the standard deviation).

Figure 3 .
Figure 3. Power spectral densities (PSD) of the bead position fluctuations.(a) Typical PSDs measured in pure water using microstructures with designed bead radii of 2 µm, 4 µm, and 6 µm.The shaded area indicates the 20-90 Hz region used for fitting the spectra.The fitted curves are plotted as solid lines.The dashed curve shows a 1/f 2 dependence.(b) PSDs obtained in different water-glycerol mixtures using an R d = 4 µm bead.The legend indicates the solution viscosities.

Figure 4 .
Figure 4. Viscosity measurements.The hydrodynamic resistance values obtained with microstructures immersed in different water-glycerol mixtures plotted against the solution viscosity.All the measurements were repeated with three different microstructures, and the error bars represent the errors of their mean γ values.